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ABSTRACT 


Doppler tracking data of three orbiting spacecraft have been reanalyzed to develop a new gravitational 
field model for the planet Mars, GMM-1 (Qoddard }£ais Model-1). This model employs nearly all available data, 
consisting of approximately 1100 days of S-band tracking data collected by NASA's Deep Space Network from 
the Mariner 9, and Viking 1 and Viking 2 spacecraft, in seven different orbits, between 1971 and 1979. GMM-1 
is complete to spherical harmonic degree and order 50, which corresponds to a half wavelength spatial resolution 
of 200-300 km where the data permit GMM-1 represents satellite orbits with considerably better accuracy than 
previous Mars gravity models and shows greater resolution of identifiable geological structures. The notable 
improvement in GMM-1 over previous models is a consequence of several factors: improved computational 
capabilities, the use of optimum weighting and least squares collocation solution techniques which stabilized the 
behavior of the solution at high degree and order, and the use of longer satellite arcs than employed in previous 
solutions that were made possible by improved force and measurement models. The inclusion of X-band tracking 
data from the 379-km altitude, near-polar orbiting Mars Observer spacecraft should provide a significant 
improvement over GMM-1, particularly at high latitudes where current data poorly resolves the gravitational 
signature of the planet 
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1. INTRODUCTION 

Knowledge of the gravitational field, in combination with surface topography, provides one of the 
principal means of inferring the internal structure of a planetary body. By removing the gravitational signal of 
the topography, the distribution of internal density anomalies associated with thermal or compositional 
differences can be estimated. Gravity can also be used to understand the mechanisms of compensation of surface 
topography, providing information on the mechanical properties and state of stress of the lithosphere. 

The earliest global gravitational field models of Mars were derived from Doppler tracking data of the 
Mariner 9 spacecraft [Lorell et aL, 1972; 1973; Bom, 1974; Jordan and Lorell, 1975; Reasenberg et aL, 1975; 
Sjogren et aL, 1975], These models provided estimates of low degree spherical harmonic gravity coefficients that 
yielded information on the oblateness and rotation vector orientation of Mars. Later models that incorporated data 
from Mariner 9 and the Viking 1 and 2 orbiters [Gapcynski et al., 1977; Reasenberg, 1977; Christensen and 
Balmino, 1979; Christensen and Williams, 1979] resolved higher degree gravity coefficients, showing the higher 
power in the Martian gravity field compared to Earth’s, and the strong correlation of long wavelength gravity 
with topography. The subsequent inclusion of additional Doppler data by Balmino et aL [1982] resulted in what 
was then the highest resolution Martian gravitational model to date: an 18* degree and order field with half 
wavelength resolution of approximately 600 km. That field, which is characterized by a spatial resolution 
comparable to what was then the highest resolution (16x16) topographic model [Bills and Ferrari , 1978], was 
utilized in analyses of the state of stress of the Martian lithosphere and the isostatic compensation of surface 
topography [Sleep and Phillips, 1979, 1985; Banerdx et aL, 1982, 1992; Willemann and Turcotte, 1982; Esposito 
et aL, 1992], However, the resolution and quality of the current gravity and topographic fields (particulaiy the 
latter) are such that the origin and evolution of even the most prominent physiographic features on Mars, nam ely 
the hemispheric dichotomy and Tharsis rise, are not well understood. 

The resolution of the Balmino et aL [1982] gravity field was limited not by data density, but rather by 
the computational resources available at the time. Because this restriction is no longer a limitation, we have re- 
analyzed the Viking and Mariner data sets and have derived a new gravitational field, designated GMM-1 
(Goddard Mars Model-1). The objectives of this study were: (1) to develop the best possible a priori 
gravitational field model for orbit determination of the Mars Observer spacecraft in support of the Radio Science 
and Mars Observer Laser Altimeter investigations, (2) to validate analysis techniques to be implemented in Mars 
Observer gravity modeling studies, and (3) to improve scientific interpretations of geophysical and geological 
data collected in previous missions to Mars. 

GMM-1 is complete to spherical harmonic degree and order 50. The corresponding half wavelength 
resolution, which occurs where the coefficients attain 100% error, is 200-300 km where the data permit. In 
contrast to previous models, GMM-1 was solved to as high degree and order as necessary to nearly exhaust the 
attenuated gravitational signal contained in the tracking data. This was possible mainly due to the use of 
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optimum weighting and least squares collocation solution techniques [Lerch et aL, 1979], which stabilized the 
behavior of the solution at high degree and order where correlation and data sensitivities become problematic. As 
discussed later, the extension of the model to high degree and order significantly reduced errors resulting from 
spectral leakage coming from the omitted portion of the gravitational field beyond the limits of the recovered 
model. GMM-1 has a higher spatial resolution than preliminary versions of this model [Smith et aL, 1990a; 

Zuber et aL, 1991], and in addition is fully calibrated to give a realistic error estimate from the solution 
covariance. 

In the following sections we discuss the development of GMM-1 and make a detailed comparison of the 
field with the previous model of Balmino et aL [1982]. We also include an error analysis and a discussion of the 
Implications of GMM-1 for Martian geophysics and for navigation and precision orbit determination in support 
of the upcoming Mam Observer Mission. 

1.1 General Approach of Gravitational Field Recovery 

Figure 1 shows a flow chart of the procedure used in the recovery of the gravitational field model. The 
data were processed using the GEODYN/SOLVE orbit determination/estimation programs [Putney, 1977]. These 
programs have previously been used in the derivation of a series of Goddard Earth gravitational Models, GEM, 
[e.g. Lerch et al., 1979; Marsh et aL, 1988; 1990], and have been adapted for the analysis of planetary tracking 
data [Smith et aL, 1990b; Nerem et al., 1993]. GEODYN provides orbit determination and geodetic parameter 
estimation capabilities, and numerically integrates the spacecraft Cartesian state and the force model partial 
derivatives employing a high-order Cowell predictor-corrector method. The force modeling includes a spherical 
harmonic representation of the gravity field, as well as a point mass representation for the Sun, Earth, Moon, and 
the other planets. Atmospheric drag, solar radiation pressure, measurement and timing biases, and tracking station 
coordinates can also be estimated. The least squares normal equations formed within GEODYN may be output to 
a file for inclusion in error analyses and parameter estimations. The SOLVE program then selectively combines 
the normal equations formed by GEODYN to generate solutions for the gravity field and other model parameters. 
The resulting gravitational model may be input back into GEODYN for residual analyses. 

1.2 Representation of The Gravitational Potential 

The gravitational potential at spacecraft altitude, is represented in spherical harmonic form as 
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where r is the radial distance from the center of mass of Mats to the spacecraft, <j> and X are the areocentric 
latitude and longitude of the spacecraft, r M is the mean radius of the reference ellipsoid of Mars, CM is the 
product of the universal constant of gravitation and the mass of Mats, P to are the normalized associated 
Legendre functions of degree / and order m, and are the normalized spherical harmonic coefficients 
which were estimated from the tracking observations to determine the gravitational model, and N is the 
maximum degree representing the size (or resolution) of the field. The gravitational force due to Mara which acts 
on the spacecraft corresponds to the gradient of the potential, V M . 

2. REFERENCE SYSTEMS FOR DSN TRACKING OF MARS SPACECRAFT 

Spacecraft orbiting Mara are tracked from Earth through the NASA DSN (Deep Space Network) 
tracking stations at Goldstone (California), Canberra (Australia) and Madrid (Spain). The adopted planetary 
ephemeris for Earth, Mara, and the other planets was the JPL DE-96 system [Starufish et al, 1976]. 

The inertial coordinate system for the Earth is defined by the direction of the Earth's rotation axis and 
the location of the vernal equinox. The orientation of the Mara rotation axis is specified by the right ascension 
and declination of the Martian pole, as given in Davies et aL [1986; 1992]. The z-axis of the Mara inertial 
coordinate system is the instantaneous Mara rotation axis. The direction of the x-axis (the IAU vector) is defined 
to be the intersection of the instantaneous Mara equator with the mean Earth equator of the appropriate 
ephemeris epoch. For this analysis, 1950.0 was chosen as the epoch, and thus the IAU vector was defined with 
reference to the Earth Mean Equator and Equinox of 1950.0 (EME50). The prime meridian of Mara is defined in 
Davies et aL [1989]. At the beginning of this analysis, tests were performed using a reference date and planetary 
ephemeris of 2000.0 but results showed that the 1971 Mariner 9 data and 1976-78 Viking data were better 
satisfied with a reference date of 1950.0. 

Various reference system constants that were used are given in Table 1. 

3. DATA SUMMARY AND ORBITAL SENSITIVITY OF GRAVITY 
3.1 Satellite Orbital Characteristics 

The Mariner 9 (M9), Viking Orbiter 1 (VOl) and Viking Orbiter 2 (V02) spacecraft were in highly 
eccentric orbits with periods of approximately one day for Viking 1 and 2 and one-half day for Mariner 9. The 
satellite orbit characteristics are summarized in Table 2. The orbital periods are nearly commensurate with the 
rotational period of Mare (24.623 hr) which produce dominant resonant perturbations [Kaula, 1966] for all orders 
m of the Viking spacecraft (24 hour period) and for the even orders of M9 (12 hour period). The resonant 
periods range mostly from about 1 to 50 days for shallow resonant terms and also include deep (very long) 
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resonant periods. The beat period, or fundamental resonant period, identifies the shift (or "walk") in successive 
ground tracks and is useful in mapping the orbital coverage over Mars (a phis sign represents an eastward "walk" 
and a negative sign for a westward "walk"). The beat period changes after cadi maneuver of the Viking 1 and 2 
spacecraft [Snyder, 1979]. For example, an orbit maneuver by V02 on March 2, 1977 produced a very slow 
walk to synchronize with the Viking Lander (VL2). Subsequently on April 18, 1977, another maneuver resulted 
in a walk around the planet in 13 revolutions, producing a beat period of 12 days. Strong resonant perturbations 
of long period were produced on VOl commencing on Dec. 2, 1978 to provide a slow walk around the planet 
Mission events such as leakages, attitude control jetting, and other phenomena that cause variations in the orbital 
periods are described by Snyder [1977; 1979]. 

The 300-km periapsis altitude orbits of VOl and VQ2 provide the strongest contribution of data to the 
solution for the higher degree terms, particularly the VOl low orbit with a range of 180° for the argument of 
periapsis (on) as compared to 25° for the V02 low orbit The observing period for the VOl low orbit shown in 
Table 2 covers almost 2 years from 77-03-12 to 79-01-27, and the periapsis point varies in latitude from +39° to 
-39° during this period. The VOl low orbit provides about a 9° ground track "walk" per revolution for the 1 1/2 
year period from 77-07-01 to 78-12-02. This corresponds to a "near repeat" of the ground track for a 39-day 
period. After the-39 day near repeat period the orbital ground track shifts by 1.8° from the previous repeat track 
which corresponds to a deep orbital resonance with a period of about 200 days. This produces, over a 200-day 
coverage, a global grid ( 39 0 latitude) with approximately 1.8° ground track separations and provides for a high 
resolution recovery of the gravity field. 

3.2 S-band Doppler Tracking Data Used in the Solution 

3.2.1 Data Summary 

The data set consisted of 263 orbital arcs representing over 1100 days of S-band Doppler tracking data 
from the Mariner 9 and Viking 1 and 2 spacecraft, collected by the Deep Space Network between 1971-1978. 
These data, grouped by satellite periapsis altitude and inclination, are summarized in Table 3. In total over 
230,000 total observations were included in the GMM-1 solution. 

3.2.2 Data Characteristics 

The data consist of two-way S-band (22 GHz) Doppler measurements compressed to 60 seconds (1 
minute data points). Data far removed from periapsis (approximately greater than 12,000 km altitude) were 
compressed over 10 minute intervals. 

All observations were collected by three DSN sites located at Goldstone, Madrid, and Canberra. They 
were processed in the differenced-rangc Doppler formulation taking into account relativistic bending due to the 
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Sun [Moyer, 1971]. Observations near satellite periapsis are most valuable for determining the gravity field and 
periapsis is generally observable by at least one of the DSN sites except when occulted by Mars. The data 
distribution and coverage per satellite orbit is discussed in Section 3.4. 

The signal is significantly degraded in precision during solar conjunction due to the solar plasma effects 
when its path comes within 5° of the Sun from Mars. This occurred for a period of about 1 months beginning 
November 7, 1976 for seven Viking 1 and 2 arcs. The data were downweighted in the solution for this period. 

33 Spectral Sensitivity of Gravity Signal 

The spectral sensitivity of the gravity field is analyzed for the Mariner 9, VOl and V02 orbits. 
Sensitivity for the high degree terms (>30) is the main area of interest and these arc compared with a threshold 
level corresponding to the precision of the DSN signal. The signal when compressed to 1 minute data points has 
a precision of 1 mm s* 1 and approximately 0.3 mm s' ! for 10 minute data points. A sensitivity study for the 
above Mars orbiters has been made by Rosborough and Lemoine [1991] and Lemoine [1992] for terms through 
degree 20. Analysis for the high degree terms is discussed in detail in Lerch et al. [1993]. A brief s ummary is 
given here. 

The orbit perturbations were studied using linear perturbation theory, and through numerical intergration 
by GEODYN. The gravity signal for sensitivity analysis employed a form of Kaula’s rule, 13x10 ‘ 5 /P f for terms 

of degree /, which was obtrained by Balmino et al [1982] for the power spectrum of Mars. The velocity 
perturbations were compared with the noise of the DSN Doppler data. Both the analytical and numerical studies 
confirm the importance of the resonance perturbations in determining the satellite sensitivity to the Mars gravity 
field. 

The resonances on the Viking spacecraft fall into three classes: (1) resonances at the low orders 
(characterized by periods of up to 40 days), (2) long period resonances (periods greater than 50 days) at specific 
higher orders, and (3) intermediate resonances at the other orders. The long period resonances result from a near 
repeat of the ground trace after an integer number of spacecraft revolutions. Thus, referring to Table 2, the 
Viking 1 orbit from July 1, 1977 through December 2, 1978, the near repeat of the ground trace (to within 1.8°) 
after 38 revolutions produces a perturbation at older 38 with a period of about 200 days. 

The analytical velocity spectrum by degree is presented in Table 4 for both M9 and the 300-km, 800- 
km, and 1500- km VOl and V02 orbits. The analytical velocity spectrum is obtained by computing the Kepler 
element perturbations using Raula's [1966] theory, and then mapping these to velocity space. Since we are 
interested in the satellite sensitivity to the gravity field over the periods of the arc lengths of data used in the 
GMM-1 solution, perturbations with periods greater than 40 days were excluded. In addition, those perturbations 
with periods between eight and 40 days have been prorated to eight days by the factor 8/period. The VOl 300- 
km orbit has a sensitivity in excess of 1 mm s' 1 (the accuracy of the S-band data) out to degree 50. In contrast. 
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the VQ2 300-km orbit is sensitive only to terms out to approximately degree 30. As the periapsis altitude is 
raised, the sensitivity in degree is diminished. The limit is degree 18 for die V02 800-km orbit, and degree 11 
for the VOl 1500-km orbit The M9 orbit has stronger pertuibadons than the VOl 1500-km orbit by virtue of its 
closer average distance to Mars, with its twice per day revolution, and smaller orbital eccentricity. 

The sensitivity was also evaluated through numerical integration using the GEODYN program for the 
VOl 300-km orbit The spectral nns velocity perturbation by order is shown in Figure 2 for different arc lengths. 
The results show sensitivity greater than 1 mm s' 1 for the high degree and order terms for arc lengths greater 
than three days. For arcs of three days, the limit in sensitivity is approximately order 30, whereas for the one day 
arcs the limit in sensitivity is approximately order 20. The increase in sensitivity results from the sampling of the 
medium period resonance perturbations. Although these results suggest it would be beneficial to process the VOl 
300-km data in batches of 8 to 16 days, this was not possible because of insufficient tracking coverage and 
errors in the nonconservative force models. 

For the highly eccentric Viking cubits, the sensitivity of spherical harmonic coefficients depends not 
only cm the periapsis altitude, but also on the location of the argument of periapsis. The GEODYN spectrum of 
rms velocity perturbations sampled by degree and order are given for an eight day arc for two VOl 300-km 
orbits in 1978. In the first arc, beginning January 13, 1978, periapsis is located near the equator (©~175°). In the 
second test arc, beginning December 20, 1978, the periapsis is located near 39°S (©-269°). When © is near 180°, 
the orbit tends to be sensitive to terms of high degree and high order, whereas when © is near 270°, the orbit is 
sensitive to terms of high degree and low order (see Figure 3). 

Another important characteristic of these eccentric orbits is that significant sensitivity to the high degree 
terms exists over a broad range of altitudes. As a demonstration, for the VOl arc described in Figures 2 and 3 
(epoch 01-15-78), the perturbations for terms of order 25 (degrees 31 to 50) are shown over a revolution in 
Figure 4. For these terms, significant sensitivity is apparent up to an altitude of 10,000 km, covering half of an 
orbit revolution. 

hi summary, the high degree sensitivity of the Viking Orbiter tracking data to the gravity field of Mars 
is determined by the periapsis altitude, location of the argument of periapsis, and the length of the arcs used to 
process the tracking data. 

3A Distribution of Observational Coverage 

The groundtrack for each of the Doppler observations is plotted in Figures 5a and 5b for a complete set 
of ground tracks covering all major data sets used in the solution. The separation between ground tracks for the 
orbital data sets is indicated in the figures by the term "walk". Also the data span is given along with the "walk" 
to depict the extent of the coverage over all data of this type as originally given in Table 2. hi these figures we 
can see the extent of periapsis coverage of the VOl low orbit ( 39 0 latitude) which is well complemented in the 
northern hemisphere by the V02 low orbit. These figures show reasonably good global data coverage for the 
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VOl low orbit, V02 low orbit, V02 800 km orbit, and also for M9. The global coverage provides for good 
separability of the lower degree terms of the gravity field and possibly out through degree 30 considering the 
strong sensitivities to these terms. Figure 6 shows the combined coverage of the low orbits of VOl and V02 
from observations with altitudes less than 5000 km and it shows coverage by different levels of altitude over 
300-km. This low altitude data coverage of the observation points along the ground tracks for the VOl low orbit 
is seen to be complemented by the V02 low orbit, particularly in the Northern Hemisphere. However, the lack of 
complete data coverage near periapsis indicates that separability will not be complete for the high degree terms 
(30 to 50) as noted above in the argument of periapsis coverage for V02. 

Nevertheless, the result from Figure 6 indicates that there is great sensitivity to the higher degree terms 
for altitudes less than 2000 km for the low altitude Viking orbits. Hence we may expect that the ground track 
coverage for the combined VOl and V02 low orbits, particularly for the observed coverage with altitude less 
than 500 km over a wide area, will provide for good resolution of localized geophysical features in the vicinity 
of these ground tracks. 

4. MODELING 

4.1 Physical Model 

Because the Viking and Mariner data do not provide uniform spatial coverage of Mars, the application 
of a priori constraints was critical to the development of a high degree and order solution. The Viking and 
Mariner data were initially processed using the gravity model of Balmino ex aL [1982]. However, in the final 
iteration to produce GMM-1, an intermediate solution, MGM-635, was used as *hr q priori model. 

The gravitational effects of the Martian solid body tide were included in the satellite force model, and a 
value of kj -0.05 was adopted [Christensen and Balmino, 1979]. The effects of atmospheric drag were 
incorporated into the satellite force model using a spherical model for the satellite body and the atmospheric 
density model developed by Culp and Stewart [1984]. The coefficient of drag, C^, was adjusted once per data 
arc, except for the VOl and V02 low periapsis orbital arcs, where C D was adjusted once per day. Solar radiation 
forces were calculated using a spherical model for the spacecraft body, and adjusting a reflectivity coefficient, 

C„ once per data arc. 

The solar flux at Mars at a given time was scaled from the Earth value for that date to the actual 
distance of the spacecraft from the Sun. One range rate bias was estimated few each tracking station per arc. The 
measurements were corrected for tropospheric effects using die Hopfield model [Hojfield, 1971]. The tracking 
data records did not contain meteorological data or tropospheric corrections, thus die corrections were computed 
a ss u min g standard temperature, humidity, and pressure, scaled to reflect the station height above sea level. 
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Third body gravitational perturbations on the spacecraft were computed from the point mass 
gravitational forces due to the Sun, the Earth-Moon system, the other planets, and Phobos, one of Mars’ natural 
satellites. In addition, Geodyn was modified to read an cphemeris for Phobos and to add the point mass 
gravitational acceleration due to Phobos to the total spacecraft acceleration. The ephemeris of Phobos was 
prepared at GSFC by processing optical measurements obtained from the Mariner 9 and Viking Orbiter images 
of Phobos [Duxbury and Callahan, 1988; 1989], 


4.2 Method of Solution 


4.2.1 Least Squares with A Priori Constraints 

The method of solution is a modified least squares process [Lerch et aL, 1979; Schwartz, 1976; 1978] 
which minimizes the sum (Q) of signal and noise as follows 


<? = e-^4^*ee^/. 
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(4.1) 


where the signal is given by C^, which are the normalized spherical harmonics comprising the solution 
coefficients. The parameter a,- 13x1ft 5 // 2 is the nns of the coefficients of degree / (a priori power rule) and is 
introduced to permit solutions to degree and order SO. This expression, which is based upon Kaula’s rule [Kaula, 
1966], has been obtained by Balmino et aL [1982] and represents the power in that gravity model. The noise 
given by is the observation residual (observed-computed) for the i* observation of satellite tracking data set 
type k, o k is the rms of observation residuals of data type k (generally significantly greater than the a priori data 
precision), and/ k is a downweighting factor to compensate for unmodeled error effects for each data type k 
(ideally for pure noise). 

The optimum weighting method estimates the combined weights directly, namely 
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When m i n i m i z i n g Q above using the least squares method, the normal matrix equation and error covariance is 
obtained as follows: 


N x - R 

where k is the solution, N is the normal matrix, R is the vector of residuals, and 


(4.3) 


8 


V = N' 1 , N = ^w k N k 


( 4 . 4 ) 


is the approximate form for the variance-covariance error matrix which must be calibrated by adjusting the 
weights. N k is the contribution for each satellite data set k to the normals, where it - 0 corresponds to the normal 
equations for the satellite a priori coefficient constraints for which N 0 is the matrix of Kaula constraints and the 
weight w Q is fixed at unity for the constraints. 

The process of m i n i miz i n g both signal (by application of the Kaula power rule constraints) plus noise in 
(4.1) is also known as collocation [Moritz, 1978]. The constraints bias the coefficients towards zero where they 
are poorly observed. With the conventional least squares approach (noise-only minimization) there is a problem 
of separability due to the strong correlation between many of the high degree coefficients. The absence of 
collocation (w 0 « 0 in (4.4) for GMM-1) results in excessively large power in the adjustment of the potential 
coefficients as in Figure 7. Hence, we see the benefit of the constraints which permit resolution of the high 
degree terms wherever the data permits and provide control of aliasing in the solution. 


42J2 Data Weighting and Error Calibration 

The weighting technique and error calibration [Lerch ex ai t 1988; Lerch, 1991] of the solution 
(equations 4. 1-4.4) is based upon subset solutions. The subset solution (CJ is formed by deleting a major data 
set k from the complete solution (C). The weight w k is adjusted as in equation (4.7) below by requiring that 


IACJ = K k a (AC*) 


( 4 . 5 ) 


where is an error calibration factor which ideally should equal to unity, 

IAC,| - {£(C - C.jf' 

° ( AC .) * {E O' - <*P 


( 4 . 6 ) 


and where a k 2 and o 2 are respectively the variance of the subset and the complete solutions. The sum in (4.6) is 
over all the coefficients, and the scale factor is needed for the errors since die error covariance in (4.4) is 
only an approximation [Lerch, 1991]. 

The new weights, w k ', should be adjusted so that each K k converges to 1 for all i fc, and the new weights 
are computed from 
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(4.7) 
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The process is iterated by forming a new complete solution and subset solutions from the new weights, and this 
process may continue until the weights converge. 

In a case where two solutions are based upon independent data, then (in the above notation) for a single 
coefficient parameter the two estimates give 


E(C - C) 2 = a 2 + a 2 


(4.8) 


whereas in Table 4.1 the data for the subset solution is wholly embedded in the complete solution in which case 

E(C - Cf = a 2 - a 2 ( 4 -9) 

as indicated by (4.6). This means that in our case the covariance between the square of the difference of the two 
estimates of the coefficients is equal to the difference of the variances of the subset and complete solutions. 

Thus, (4.8) and (4.9) represent extremes in estimation, complete independence and complete dependence. 


5. RESULTS OF THE GMM-1 SOLUTION 


5.1 Description of Solution 

GMM-1 is a 50 x 50 spherical harmonic gravity model. There are a total of 5250 estimated parameters: 
2597 gravity coefficients plus GM, and the arc parameters. The GMM-1 gravity coefficients through degree and 
order 50 are shown in Appendix A. Calibrated accuracy estimates of the coefficients have also been obtained for 
the model. 


52 Gravity Model Tests Using Orbital Observation Residuals 

Orbital arcs have been selected from the 7 major data sets summarized in Table 3 and used to test the 
orbital accuracy of the model by fitting the DSN Doppler data. Observation residuals have been computed from 
our 50x50 model and compared with the prior best available 18x18 gravity model of Balmino et aL [1982]. 
Table 5 is a compilation of orbit tests for 14 arcs. Each arc is fit using the Balmino et aL field and GMM-1. For 
all 14 test arcs, the RMS residual fits are significantly smaller, sometimes 5 to 10 times smaller, when computed 
using GMM-1 than when using the Balmino et aL field Table 6 summarizes the results of orbit prediction tests 
for a subset of the arcs in Table 5. The orbits obtained in the orbit accuracy tests are projected forward in time 2 
or 3 days. Then RMS residual fits are compared for the data in the predicted time periods. Again, the fits are 
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significantly smaller when computed using GMM-1 than when using the Balmino et al. field The improvements 
in the fits is not entirely due to the increased resolution of GMM-1. An 18x18 version of GMM-1 outperformed 
the 18x18 model of Balmino et aL in all cases except the 300-km VOl and V02 orbits, for which the 
performance was comparable. 

S3 Analysis of the Gravity Coefficients 

Figure 8 is a plot of the degree variance of the coefficients (power) and error variance of GMM-1 per 
degree. Also plotted are the power of the 18x18 gravity field [Balmino et al. t 1982] and a power rule (13x10 s // 2 ) 
taken from Balmino et aL [1982], which is the basis of the constraint matrix used in GMM-1. The plot shows 
that for degrees less than 15, the power spectrum of GMM-1 and the Balmino et aL field are about the same. 
However, above degree 15, GMM-1 drops below while the Balmino et aL field rises above the power spectrum 
of Balmino’s rule. The upward turn of the Balmino et aL field is undoubtedly due to aliasing. Aliasing adversely 
affects the performance of a gravity field with respect to orbit fits from independent data, orbit predictions, and 
other geophysical information which are derived from the gravity coefficients. The truncation level of GMM-1 at 
degree 50 is high enough so that the high degree gravity signal is not significantly aliased into the lower degree 
terms. 

The power spectrum of GMM-1 drops below the values of the power rule for high degrees. Above 
degree 22, the errors of the coefficients are larger than the coefficients themselves. This drop off in the power 
spectrum occurs because the drag parameters (once per day values) are absorbing part of the gravity signal. Also, 
the high degree terms are highly correlated and hence the effect of the power rule constraint in the solution is 
quite strong which further explains the small power for these terms. However, because the a priori constraint 
does not have a major effect on the solution, the terms do contain information on the short wavelength gravity 
field in the vicinity of the spacecraft periapses. While the power in the field falls below that predicted by the 
power rule at high degrees, the field is a better representation of the true gravitational signature of the planet at 
those wavelengths than would be the case if the field were solved to lower degree and order and all of the high 
degree and order coefficients were constrained to zero. 

Figure 9 shows the coefficient differences between GMM-1 and the 18x18 gravity field. While the rms 
differences per degree are about the size of Balmino's rule for terms above degree 10, the differences between 
particular ( C^, pairs even for lower degree terms are seen to be quite large. In fact, the rms differences for 
lower degree terms are over an order of magnitude greater than the error estimates of GMM-1 as given in Figure 
8. The coefficient discrepancy between these models reflects the large differences seen in the orbital residuals for 
the two models as shown in Table 5. 

SA Calibration of Gravity Model Errors 
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The calibration of the gravity model error estimates is based upon the method described in Section 422 
a pd is developed in greater detail by Lerch et aL [1991]. In the application of this method, weights of basic 
observation sets from different orbits are adjusted based on subset solutions. The data is separated into 7 groups 
(see Table 6) yielding 7 subset solutions for the weight adjustment In Table 6 each group is assigned an a priori 
data weight which is based on our experience in computing previous gravity solutions. For example, the Viking 1 
1 50Q-fan data group is assigned an error of 1 cm s' 1 while Viking 1 300-km data group is assigned an error of 
.71 cm s 1 [wt-l/(.71)2-2]. The larger errors (indicating down-weighting) for the data sets of Viking 1 at 1500- 
km and Viking 2 at 1500-km (55° inclination) are due to the synchronous (repetitive) nature of the orbits (over 
the V ikin g Landers) as shown in Table 2 and in Figure 5a and 5b for the data distribution. The calibration 
factors (Jfc) given in Table 7 indicate that the model is reasonably well calibrated where a factor of *-l indicates 
perfect caUbration. 

5.5 Error Analysis 

The error covariance matrix, which is calibrated in Section 5.4, was used to project the orbital errors in 
satellite position and velocity. Table 8 shows the projected errors for M9, V01-300 km, VO @-300 km, and MO 
for a 6-day arc length. The results for the Mars Observer orbit are of special interest since these errors will affect 
the orbit determination. Figures 10 and 1 1 give respectively the error spectrum by degree and order for the radial 
and along -track position components of Mars Observer (cross-track errors are similar to the radial errors). Note 
the largest error is shown for resonant order 25 indicating that a field complete to at least degree 30 is required 
to reasonably model these coefficients based upon the error spectrum by degree. 

5.6 Recovery of GM 

In the GMM-1 solution, the GM of Mars was adjusted along with the other coefficients of the Mars 
gravity field. The value of GM determined in the solution was 42828.36 0.05 km V 2 . Lemoine [1992] analyzed a 
smaller set of Viking and Mariner 9 Doppler data as well as Viking Orbiter range data and determined a value of 
CM of 42828.40 0.03 km 3 s‘ 2 . The estimates of the Mars GM from Mariners 4, 6 and 9 [Null, 1969; Anderson et 
aL f 1970; O'Neill et a/., 1973] are in close agreement with the GMM-1 value. The Mariner 4, 6 and 9 values are 
especially interesting since they are derived from tracking of spacecraft from a flyby of the planet Mars. In these 
cases, the estimate of the Mars GM is largely uncoupled from the remaining coefficients of the Mars gravity 
field. 

6. GEOPHYSICAL IMPLICATIONS 

Figure 12 shows free air gravity anomalies computed from GMM-1 complete to degree and order 50, 
and Figure 13 displays accompanying gravity anomaly errors computed from the error covariance matrix. As 
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illustrated in Figures 5 and 6, the satellites used in this study are characterized by a complicated distribution of 
low altitude data. The shortest wavelengths resolved (200-300-km half wavelength) occur within the latitudinal 
band of 40 0 corresponding to the data coverage from the VOl low orbit This region also includes the periapsis 
coverage (0° to 30° latitude) of the VQ2 low orbit as seen in Figure 6. That figure also shows that above 40° N 
latitude there is still strong coverage of periapsis extending to 63°N latitude particularly for the VQ2 800-km 
altitude orbit This coverage is reflected in the gravity anomaly error map which shows more longitudinal 
structure and better resolution than in the corresponding southern hemisphere beyond the region of 40°S latitude. 

In Figure 12, the free air anomalies are overlain by contours of topography. The topographic field is a 
spherical harmonic expansion, also complete to degree and order 50, of the Mars Digital Elevation Model [DEM; 
Wu et aL, 1986]. The spherical harmonic topographic model was defined to have zero mean elevation, and so 
while the spherical harmonic and DEM have similar hypsometric distributions, elevation values from the former 
are offset by approximately two km from the latter. 

As in previous studies, the gravity anomalies correlate well with principal features of Martian 
topography, including volcanic shields, impact basins and the Valles Marineris. Most major features exhibit 
anomalies with considerably higher magnitudes than in previous models. GMM-1 also exhibits gravity anomalies 
in association with some observed structures that were not previously detected. For example, GMM-1 resolves all 
three Tharsis Montes, while the model of Balmino et aL [1982] fails to resolve the central volcano in the line, 
Pavonis Mons. GMM-1 also shows considerably more detail associated with the Valles Marineris, and for several 
of the major impact basins including Isidis and Argyre. However, it is important to interpret short wavelength 
features resolved in the model with caution, as the the coefficients associated with the highest degree and coder 
terms are 100% in error. 

One of the most prominent physiographic features cm Mars is the hemispheric dichotomy, which is 
characterized by a 2 -km elevation difference between the northern and southern hemispheres of Mars. However, 
the dichotomy does not have a distinct gravitational signature associated with it This indicates that the 
dichotomy boundary is isostatically compensated at the resolvable wavelengths of GMM-1, perhaps due to a 
change in crustal thickness across the boundary, as suggested in previous studies [Phillips and Saunders, 1975; 
Lambeck, 1979; Phillips and Lambeck, 1980; Phillips , 1988]. 

It is significant to note that several prominent anomalies in GMM-1 fail to correlate with observed 
surface features. These include a 300-mgal negative anomaly on the western edge of Tharsis (lon-200°E, 
lat-+20 |O N) and a 200-mgal positive anomaly in Utopia (loo- 105 °E, lat"50°N)* Both of these areas are in the 
nothem hemisphere and may have been resurfaced. These features were also present in the field of Balmino et 
aL [1982], but the anomalies were smaller in magnitude. 
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As for the gravity anomaly representation of the field, the geoid from GMM-1 as shown in Figure 14 
exhibits a higher dynamic range of power (2300 m vs. 1950 m) than the model of Balmino et aL [1982]. The 
distribution of geoid errors shows a similar pattern to the gravity anomaly errors. 

A detailed geophysical interpretation of GMM-1, which includes a spectral analysis of the gravity and 
topography fields and a global inversion of die fields for simultaneous estimations of density anomalies in the 
Martian crust and mantle, is presented in a companion paper by Bills et aL [manuscript in preparation, 1993]. 

1. SUMMARY 

Re-analysis of Doppler tracking data from the Mariner 9 and Viking 1 and 2 spacecraft has led to the 
derivation of a 50* degree and order gravitational model for Mars. The model has a maximum (half wavelength) 
spatial resolution of 300-km where the data permit, which represents a factor of two improvement over that 
attained by the previous field of Balmino et aL [1982] which utilized essentially the same data. Probable reasons 
for the significant improvement achieved include: increased computational capabilities, the application of 
collocation and optimum data weighting techniques in the least squares inversion for the field, and the use of 
longer arcs (days vs. hours) than used previously for Viking low altitude data made possible by improved fcoce 
and measurement models. 

Error analyses based on the observation data, derived power spectrum, and comparison with topography 
demonstrate that this field represents the orbits with considerably better accuracy and shows greater resolution of 
identifiable geological structures than previous models. The model also shows a greater dynamic range of power 
in both the gravity anomaly field and the geoid. The inclusion of X-band tracking data from the 379-km 
altitude, near-circular, polar orbit of Mars Observer will allow significant improvement of the Martian 
gravitational field [Smith et aL, 1990b; Esposito et aL,\99Q\ Tyler et aL, 1992], with the greatest refinement 
occurring at high latitudes far removed from the Mariner 9 and Viking 1 and 2 periapsis latitudes. That 
gravitational field, in combination with topography data from the Mars Observer Laser Altimeter (MOLA) 

[Zuber et aL, 1992] will allow detailed analyses of Mars’ internal structure, state of lithospheric stress, and 
mechanisms of isostatic compensation of surface topography. 
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Table 1. Planetary and Astronomical Constants 


Parameter 

Value 

Unit 

Gravitational Constant (GM) 

42828.28 

km 3 s' 2 

Equatorial Radius of Mars 

3394.2 

km 

Spin Rate of Mars 

350.891983 

deg day* 1 

Solid Tide Amplitude (k^ 

0.05 


Flattening of Mars 

1/191.1372 


Speed of Light 

2.99792458x10 s 

m s' 1 

Astronomical Unit 

1.49597870660xl0 u m 
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Table 2. SATELLITE ORBIT CHARACTERISTICS 
Including Beat Periods and Ground Track Walks 


Comments 

New orbit U 

End of data I 

Synchronous 1 

9 rev. near repeat 
Over lander 
Near Phobos 
New orbit 
New walk 
Dual station 
Near synch. 

End of data 

New orbit 
Synchronous 
New Ind. 

Low perlaps.And. 
change 
Over lander 
13 rev. near repeat 
Near Detmos 
Near Dei mas 
Low perlaps. 

End of data 

Walk per 

revolution 

(degrees) 

% 

o 

r-* 

<1 

44 

<1 

25 

43 

17 

10 •• 
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8 
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27.31 
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26.79 

26.48 

24.73 

22.72 

24.29 

24.20 

23.98 

Inclination 

(degrees). 
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64. 0.62 
64, 0.60 
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Walk for two 12 hr. revolution for Mariner 9 orbits 
Start epoch of the orbit parameters cited 

Secondary walk is about 2 degrees over approximately 200 day duration 



Table 3. 

Summary of Data used In Mars Gravity Model GMM-1 
DSN Tracking Data (+ .1 cm / see) 
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Tabic 4, 


HAP Velocity Perturbations By Degree 


USING A POWER RULE OP 13E-0S/L"2 
MARINER 9, VIXING 14 2 SAMPLED ORBITS "* 
VELOCITY PERTURBATIONS IN CH/SEC 


FERIAPSIS 

vxci 

VXG2 

VXG2 

VXG1 

KRN9 

ALTITUDE (XH) 

i 300 

300 

• 00 

• 1500 

1500 

DEGREE 

EPOCH 

EPOCH 

EPOCH 

EPOCH 

EPOCH 


78-01-15 

77-12-17 

77-04-19 

77-02-05 

72-04-10 

2 

544.175 

453.411 

144.452 

44.324 

828.121 

3 

353.495 

218.135 

39.741 

34.417 

100.049 

4 

207.188 

104.429 

14.519 

13.981 

55.852 

S 

107.358 

59.410 

7.771 

5.459 

24.812 

4 

44.205 

34.485 

4.248 

2.270 

12.443 

7 

40.783 

21.119 

2. 498 

1.051 

4.003 

8 

28.410 

13.923 

1.417 

0.579 

2.449 

9 

18.054 

9.347 

1-397 

0.304 

1.135 

10 

13.174 

4.477 

1.118 

0.158 

0.490 

11 

9.049 

4.734 

1.523 

0.093 

0.240 

12 

4.108 

3.555 

2.484 

0*054 

0.147 

13 

4.711 

2.449 

1.195 

0.031 

0.105 

14 

3.858 

2.053 

1.350 

0.019 

0.041 

18 

1.390 

0.780 

0.211 

0.010 

0.002 

22 

0.581 

0.3B1 

0.048 

0.002 

0.000 

24 

0.317 

0.249 

0.012 

0.000 

0.000 

30 

0.240 

0.110 

0.004 

0.000 

0.000 

34 

0.370 

0.047 

0.000 

0.000 

0.000 

34 

0.507 

0.034 

0.000 

0.000 

0.000 

38 

0.448 

0.028 

0.000 

0.000 

0.000 

40 

0.373 

0.029 

0.000 

0.000 

0.000 

42 

0.241 

0.034 

0.000 

0.000 

0.000 

44 

0.245 

0.038 

0.000 

0.000 

0.000 

50 

0.322 

0.027 

0.000 

0.000 

0.000 


* HAP; HARMONIC ANALYSIS OP PERTURBATIONS PROM ANALYTIC THEORY. 

* PERIODS GT. S DAYS HAVE THEIR AMPLITUDES MULTIPLIED BY { 8/PERIOD). 
PERIODS CT. 40 DAYS HAVE BEEN EXCLUDED. 


Table 5, Orbit Accuracy Tests : RMS of Orbital Fits in cm/sec 


Arc 

« 

Satell Ite 

Arc 

epoch 

yymmdd 

No. of 
obs. 

Arc 

length 

days 

Balmlno 
18 x 18 

GMM-1 
50 x 50 

1 

Mari ner-9 
Inc. = 64° 

7201 13 

1896 

B 

.456 

.090 

2 

VOl 1500 km. 
Inc. a 39° 

760822 

1326 

fl 

.687 

.097 

3 

V02 1500 km. 
inc. = 55° 

760917 

151 1 

fl 

.387 

.196 

4 

V02 1500 km. 
Inc. * 75° 

761026 

1350 

fl 

.649 

.340 

5 

V02 800 km. 
inc. * 80° 

770102 

682 

B 

.434 

.143 


Average 

.522 

.173 


6 


771 122 

568 

9 

5.07 

1.04 

7 

VOl 300 km 

780210 

754 

9 

6.44 

1.22 

8 

Inc. = 39° 

780604 

538 

2 

2.42 

.74 

9 


78081 1 

387 

2 

1.43 

.09 

10 


780904 

1025 

8 

8.58 

1.24 


Average 

4.79 

.87 


11 


771 1 17 

11 14 

6 

1.52 

1.02 


V02 300 km 






12 

inc. x 80° 

771217 

688 

2 

.73 

.11 

13 


780516 

791 

8 

7.67 

1.78 

14 


780526 

705 

4 

.60 

.15 


Average 

2.63 

.77 


For arcs 1-5 arc parameters adjusted are position, velocity , C r and station biases 
For arcs 6- 14 arc parameters adj us ted are posi tlon , velocity , C f and Cj per arc 
Arc 13 comes in as 2 separate arcs in GMM-1 
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Table 8 . Orbit Prediction Tests : RMS of Prediction Fits In cm/sec 


Arc 

# 

Sateil i te 

Arc 

epoch 

yymmdd 

No. of 
obs. 

Predict 

ed 

peri od 
days 

Balmlno 
18 x 18 

GMM-1 
50 x 50 

1 

Mari ner-9 
inc. * 64° 

7201 13 

1308 

3 

5.61 

.36 

2 

VOl 1500 km. 
Inc. ■ 39° 

760822 

640 

3 

4.20 

2.58 

3 

V02 1500 km. 
inc. * 55° 

7609 1 7 

720 

3 

1.62 

.68 

D 

V02 1500 km. 
Inc. - 75° 

761026 

697 

3 

8.04 

2.38 

5 

V02 800 km. 
inc. » 80° 

770102 

367 

3 

21.40 

.73 


Average 

8. 17 

1.34 




VOl 300 km. 
inc. - 39° 

780904 

353 

3 

105.6 

13.2 


IKS 

V02 300 km. 

771 1 17 

597 

3 

102.9 

30. 1 

KS 

inc. - 80” 
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Table 7. Calibration and Data Weights of GMM- 1 


Subset 

Solution 

Dataset 

Removed 

Apriori 
Sigma 
Weigh ts 

y 

cm / sec 

GMM' 1 
Sigma 
Weigh ts 

y 

cm / sec 

GMM - 1 
Callbratio 
n 

Factors 
k •• 

V02 1500 km 55° Inc. 

1.0 

4.1 

1.2 

V02 1500 km 75° Inc. 

1.0 

1.5 

.81 

V02 800km 80° Inc. 

.71 

.72 

.81 

V02 300 km 80° Inc. 

.71 

1.0 

1. 16 

VOl 1500 km. 39° 
Inc. 

1.0 

3.5 

.96 

VOl 300 km 39° Inc 

.71 

.8 

1.05 

Mari ner-9 

1.0 

2.0 

.99 


k = 


_L 

2 \ Z 

s( c - c ) 
^X(cf 2 - cr 2 ) j 




k=l for complete convergence 


33 







Table 8. Projected Gravity Orbit Error from GMM-1 Covariances 
( Long period terms excl uded ) 

Orbit Position Error in meters 


Satellite Length [ day si Radial Along-Track Crogg-Traclt Total 


Mars Observer 

6 

07 

757 

90 

705 

Marl ner-9 

6 

2 

4 

2 

5 

Viking- 1,300 km 

6 

26 

69 

31 

80 

Viking-2,300 km 

6 

26 

83 

12 

87 


Orbit Velocity Error in cm/sec 

Satellite Length fdavsl Radial Along-Track Crft.gg-.Track Total 


Mars Observer 

6 

66.3 

8.0 

8.0 

67. 1 

Marl ner-9 

6 

.07 

.03 

.02 

.08 

Viking-1,300 km 

0 

1.9 

.8 

.2 

2. 1 

Viking-2,300 km 

6 

2.0 

.9 

.1 

2.3 
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TRACKING DATA 
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IAVITY FIELD ERRORS) 




• GEODYN: Numerically integrated perturbations 
HAP: Harmonic analysis of perturbations from analytical theory 

Figure 2. Spectral Sensitivity of Gravity Signal (by Order) 
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(cl>= 175°) 


Velocity Perturbations In .001 cm/see 

RSS* DEC VTC1 i 78-01-15 

392339 2 ••• Xrg . of periapsis ■ 175 deg 

124393 9 •••**• 9 day arc length 

29904 10 <»•*•••* 

9027 14 54******156 

2909 II 6 52***943 29 

1297 22 3 10139592237 9 


475 

25 

1 

2 

7151299 

76 

3 








400 

30 

1 

1 

2 

20137143 

29 

2 







258 

34 ’ 

0 

0 

1 

2 

33104 

77 

13 

1 






229 

35 

0 

0 

0 

1 

13 

69 

96 

29 

2 

1 





203 

39 

0 

0 

0 

0 

4 

39 

79 

48 

9 

2 

1 




197 

40 

0 

0 

0 

0 

1 

20 

63 

62 

22 

7 

1 

0 



179 

42 

0 

0 

0 

0 

0 

9 

43 

67 

41 

19 

5 

1 

0 


195 

45 

0 

0 

0 

0 

0 

1 

14 

49 

73 

59 

29 

9 

1 0 


199 

50 

0 

0 

0 

0 

0 

0 

3 

22 

69 

84 

67 

34 

11 0 

0 


ORDt 

2 

5 

10 

14 

19 

22 

25 

30 

34 

36 

39 

40 

42 46 

50 

* RSS taken 

over all 

orders 

t , not 

juet sampled orders. 



(ti>=269°) 


Velocity Perturbations In .001 cb/mc 


RSS* 

DEC 



\TClt 

78-12-20 




591650' 




Xrg. 

of perlapsii 

■ 1 

- 269 deg 

72711 

6 



8 

day arc 

length 




22705 

10 

*••494336 











9766 

14 


3 










4959 

19 


16 

1 









2771 

22 

762***554 

37 

6 

0 








1631 

26 

146420411 

62 

11 

1 

0 







994 

30 

120 47192 

91 

2 

2 

0 

0 






604 

34 

176108 38 

65 

13 

1 

1 

0 

0 





493 

36 

133175148 

45 

5 

2 

0 

0 

0 

0 




378 

39 

157131 40 

24 

15 

1 

0 

0 

0 

0 

0 



316 

40 

33 63 81 

46 

6 

2 

0 

0 

0 

0 

0 

0 


244 

42 

112106 60 

6 

10 

2 

0 

0 

0 

0 

0 

0 

0 

162 

46 

67 69 52 

18 

2 

2 

0 

0 

0 

0 

0 

0 

0 0 

no 

50 

32 37 35 

19 

3 

1 

0 

0 

0 

0 

0 

0 

0 0 


ORDt 

2 6 10 

14 

18 

22 26 

30 

34 

36 38 

40 

42 46 

* RSS taken over all 

orders 

, not 

just s 

anpled orders. 


Figure 3. Signal Sensitivity by Degree and Order for VOl Low Orbit 
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* based on 8-day arc from epoch 78-01-15 

Figure 4. Velocity Perturbation for VOl Low Orbit due to 25th-Order Harmonic 
Coefficients for (a): Degrees 31-50, (b) Degrees 41-50* 
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Figure 6 

Viking- 1 and Viking-2 300 km. Observations with Altitude <5000 km. 
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Figure 9. Coefficient Differences for Balmino minus GMM-1 
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RMS per Degree 


Total RMS = 757m 


RMS per Order 



Harmonic degree 



Harmonic order 


Figure 10. Projected Radial Position Error on Mars Observer from GMM-1 
Gravity Covariances 
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Figure 11. Projected Along-Track Position Error on Mars Observer 
from GMM-1 Gravity Covariances 


Total RMS * 757m 

RMS per Degree RMS per Order 



Harmonic degree Harmonic order 


Figure 11. Projected Along-Track Position Error on Mars Observer from GMM-1 
Gravity Covariances 
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Figure 13 

Mars Gravity Anomaly Errors from GMM 
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